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Abstract 

According to the real r-conjecture, the number of real roots of a 
sum of products of sparse univariate polynomials should be polynomi- 
ally bounded in the size of such an expression. It is known that this 
conjecture implies a superpolynomial lower bound on the arithmetic 
circuit complexity of the permanent. 

In this paper, we use the Wronksian determinant to give an upper 
bound on the number of real roots of sums of products of sparse poly- 
nomials of a special form. We focus on the case where the number of 
distinct sparse polynomials is small, but each polynomial may be re- 
peated several times. We also give a deterministic polynomial identity 
testing algorithm for the same class of polynomials. 

Our proof techniques are quite versatile; they can in particular 
be applied to some sparse geometric problems that do not originate 
from arithmetic circuit complexity. The paper should therefore be of 
interest to researchers from these two communities (complexity theory 
and sparse polynomial systems). 
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1 Introduction 



The complexity of the permanent polynomial 

n 

per(xn, . . . ,x nn ) = ^ IJ^i<7(i) 

TGS n 1 = 1 

is one of the central open problems in complexity theory. It is widely be- 
lieved that the permanent is not computable by arithmetic circuits of size 
polynomial in n. This problem can be viewed as an algebraic version of the 
P versus NP problem |24L |7] . 

It is known that this much coveted lower bound for the permanent would 
follow from a so-called real r- conjecture for sums of products of sparse 
univariate polynomials [15] . Those are polynomials in R[x] of the form 
Yli=i Ilj=i fij( x )> where the sparse polynomials fij have at most t mono- 
mials. According to the real r-conjecture, the number of real roots of such 
an expression should be polynomially bounded in k, m and t. The original 
r-conjecture by Shub and Smale [23j deals with integer roots of arbitrary 
(constant-free) straight-line programs. 

As a first step toward the real r-conjecture, Grenet, Koiran, Portier and 
Strozecki [10] considered the family of sums of products of powers of sparse 
polynomials. Those polynomials are of the form 

k m 
i=l j=\ 

They are best viewed as sums of products of sparse polynomials where the 
total number m of distinct sparse polynomials is "small", but each polynomial 
may be repeated several times. In particular, if one can find a (kt) ^2°^ m ^ 
upper bound on the number of real roots, then it will imply the real-r con- 
jecture in the case where the number of distinct sparse polynomials is loga- 
rithmically bounded. The upper bound on the number of real roots obtained 
in [TO] is polynomial in t, but exponential in m and doubly exponential in k. 

Bounds on the number of real zeros for systems of sparse polynomials 
were extensively studied by Khovanskh [14j in his "fewnomial theory". His 
results imply an upper bound exponential in k, m and t. In this article, we 
will give a bound of order t°( k m "> , thereby removing the double exponential 
from [10J while staying polynomial in t. Moreover, our results extend well 
to some other families of functions. In particular, they extend a result from 
Avendaho [4] on the intersection of a sparse plane curve and a line. He 
gave a linear bound on the number of roots for polynomials of the form 
Si=i x° li {ax + b)^ i where a% and /3j are integers and gave an example proving 
that his linear bound does not apply for non-integer powers. Our result gives 
a polynomial upper bound for the wider family ([I]) where the polynomials 
fij are of bounded degrees and the ctij are real exponents. 
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In addition to bounds on the number of real roots, we also give a de- 
terministic identity testing algorithm for polynomials of the form ([!]). The 
running time of our algorithm is polynomial in t, in the bit size of coeffi- 
cients and of the powers (a» and exponential in k and m. Polynomial 
Identity Testing (PIT) is a very well-studied problem. The Schwartz- Zippel 
lemma yields a randomized algorithm for PIT, but the existence of an effi- 
cient deterministic algorithm is an outstanding open problem. Connections 
between circuit lower bounds and deterministic PIT algorithms were discov- 
ered in 1980 by Heintz and Schnorr |1 1| . then more recently by Kabanets and 
Impagliazzo [13], by Aaronson and van Melkebeek pQ and by Agrawal (2]. 
Recently, many deterministic PIT algorithms have been found for several 
restricted models (see e.g. the two surveys In particular, a deter- 

ministic PIT algorithm for functions of the form of ([I]) has already been given 
in |10| . Their algorithm is polynomial in t, exponential in m but doubly ex- 
ponential in k whereas we give a new algorithm which is only exponential in 
k. 

We now present our main technical tools. Finding the roots of a product 
of polynomials is easy: it is the union of the roots of the corresponding 
polynomials. But finding the roots of a sum is difficult: for example how to 
bound the number of real roots of fg + 1 where / and g are i-sparse? It is 
an open question to decide if this bound is linear in t. Our main tool in this 
paper to tackle the sum is the Wronskian. We recall that the Wronskian 
of a family of functions fi, ■ ■ ■ ,fk is the determinant of the matrix of their 
derivatives of order up to k — 1. More formally, 



W(h,..., 




The Wronskian is useful especially for its connection to linear independence 
(more on this in the next section). Another classical and very useful tool is 
Descartes' rule of signs: 

Lemma 1 (Strong rule of signs). Let f = Ya=i a % xOLi be a polynomial such 
that a± < «2 < • • • < oli and ai are nonzero real numbers. Let N be the 
number of sign changes in the sequence (oi,...,Ot). Then the number of 
positive real roots of f is bounded by N. 

In this article, we will use a weak form of this lemma. 

Lemma 2 (Weak rule of signs). Let f = Ya=i a i x0ii be a polynomial such 
that «i < «2 < • • • < a t an d a-i are nonzero real numbers. Then the number 
of positive real roots of f is bounded by t — 1. 

In their book [20J, Polya and Szego gave a generalization of the strong 
rule of sign using the Wronskian. We show in Theorem [7] that bounding 
the number of roots of the Wronskian yields a bound on the number of 
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roots of the corresponding sum. In general, the Wronskian may seem more 
complicated than the sum of the functions, but for the families studied in 



this paper it can be factorized more easily (Theorems 12 and|13[). 

The paper is organized as follows. The main results of Section [2] are The- 
orems [7j [8] and [9] which bound the number of roots of sums as a function of 
the number of roots of the Wronskian. Then, in Section [3j we apply these 
results to particular families of polynomials. The main applications that we 
have in mind are to polynomials of the form ([I]), and to the polynomials 
studied by Avendaho. We give in Section [4] some PIT algorithms for poly- 
nomials of the form |T]) . The proof of Theorem [9] will be given in Section [5] 
And finally, we show in Section [6] that our method is optimal in a precise 
sense. Some of the proofs are postponed to the Appendix. 



2 Zeros of the Wronskian as an upper bound 

Let us recall that for a finite family of real functions fi, ■ ■ ■ , fk sufficiently 
differentiable, the Wronskian is defined by 



W(h,..., 




We will use the following properties of the Wronskian. 

Lemma 3. Let fi, ■ ■ ■ , fk be real functions k — 1 times differentiable. 
Then, W(gf 1: ...,gf k ) = g k W(h, ...,/*). 

As a corollary: 

Lemma 4. Let fi,- ■ ■ , fk be real functions k — 1 times differentiable and I 
be an interval where they do not vanish. 

Then, over L, we have W(f u ...,/*) = {fit W (($)'>•■•> (%)'\ 

These results can be found in |20] (ex. 57, 58 in Part 7). Notice that 
the Wronskian of a linearly dependent family of functions is identically zero 
(if a family is dependent then the family of the derivatives is also dependent 
with the same coefficients). But the converse is not necessarily true. Peano, 
then Bocher, found counterexamples |18[ QjJJ |5] (see [8] for a history of these 
results). However, Bocher |6] proved that this converse becomes true if the 
functions are analytic |12| . 

Lemma 5. If fx, ■ ■ ■ , fk are analytic functions, then (fi) is linearly dependent 
if and only if W(fi, . . . , f k ) = 0. 

Definition 6. For every function g and interval I, we will denote Zj(g) 
the number of distinct real roots of g over I. We just write Z(g) when the 
interval is clear from the context. 
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Throughout the paper we consider only intervals that are not reduced 
to a single point (and we allow unbounded intervals). The next Theorem is 



nearly implied by Voorhoeve's result (see Theorem 20) below. 



Theorem 7. Let k be a non zero integer. Let f be k functions k — 1 
times differ entiable in an interval I such that for all i < k, the Wronskian 
W(fi, . . . , fi) does not have any zero over I. 

If the real constants ax, . . . , a k are not all equal to 0, a\f\ + 02/2 + • • • + 
a k f k has at most k — 1 real zeros over I counted with multiplicity. 

Proof. We show this result by induction on k. If k = 1, then, fx = W(fx) 
does not have any zero. Moreover, ax is not zero. So, axfi has no zeros. 

For some k > 2, let us suppose that the property is true for all linear 
combinations of size k— 1. Denote z the number of zeros of ax fx + . . . -\-ctkfk- 
If 02 = «3 = . . . = flfc = 0, then a\ 7^ and ax fx + 02/2 + . . . + a k f k = 
a 1/1 has no zero, and the conclusion of the theorem holds true. Otherwise, 
eti + %p + . . . + ^j^- has z zeros (since fx = W(f%) does not have any zero 

by hypothesis). By application of Rolle's Theorem, 02 + . . -+ a k (^j^ 

has at least z — 1 zeros over /. 

Function fx does not have any root in I, so the functions (^j^j ' ' ' ' ' (ji 
are k — 2 times differentiable. Moreover, for all 2 < i < k, 
W {{h/ fx)' (fi/fx)') = W(fi, fi)/f[ does not have any roots over 
/. Since the coefficients 02, • • . , Ofc are not all zero, by induction hypothesis 

a 2 (£)+• • -+a k (ft)' has at most k-2 zeros. Hence a 1 / 1 +a 2 / 2 + . . .+a k f k 
has at most k — 1 zeros by Rolle's theorem. □ 

The following theorem gives us a method to find upper bounds on the 
number of roots. We will show in Section [6] that it is sometimes tight. 

Theorem 8. Let fx, ■ ■ ■ , fk be analytic functions on an interval I. If the 
real constants ai, . . . ,a k are not all equal to 0, 



Z(axfx + ■■■ + a k f k ) <(l + J2z(W(h, . . . , fr))) k - 1. 



(2) 



More precisely, if T = {x G I\3i < k, W(f\, . . . , fi)(x) = 0} is finite, then 

Z(ai/i + ... + a fc / fc )<(l + |T|)fc-l. (3) 

Moreover, the inequalities still hold if on the left side, zeros which are not 
zero of one of the Wronskians W(fx, ■ ■ ■ , fi) are counted with multiplicity. 

Proof. If ai/i + . . . + a k f k is the zero polynomial, then the family is linearly 
dependent and so the Wronskian W(fx, • • • , fk) is also the zero polynomial. 
This means that T = / is infinite and the inequality is verified. 
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Otherwise, a%fx + . . . + a k fk has a finite number of zeros. We have T = 

k k 

[}Z{W{h,..-Ji)). So, |T| < E\Z(W(f 1 ,...J i ))\ and we will prove (3). 

i=l i=l 

The set I\ T is an union of |T| + 1 intervals. Let J be one of these intervals. 
With Theorem[7j we get Z J (a 1 f i + . . . + a k f k ) < k - 1. So ai/i + . . . + a fc / fc 
has at most (1 + |T|) [k— 1) zeros over I\T and at most (1 + |T|) (k— 1) + |T| 
zeros over /. □ 

In Section [5] we will prove the following variation on Theorem [8] 

Theorem 9. Let fi, ■ ■ ■ , fk be analytic linearly independent functions on an 
interval I. Then, 

k-2 

Z(fi + ... + f k )<k-l + Z{W k ) + Z(W fc _i) + 2 Y, Z(Wj). 

i=i 

In most applications, this result yields a better bound than Theorems [7] 
andEl 



3 Applications 



In this section, we prove Theorem 12 which bounds the number of zeros of the 
polynomials of the form Q. The given bound improves both Grenet, Koiran, 
Portier and Strozecki's result [10] and the bound implied by Khovanskh's 
fewnomial theory [14j. At the end of this section, we also extend Avendaho's 
result to real exponents. We saw before (in Section [2J that the number of 
zeros of a linear combination of real functions can be bounded as a function 
of the number of zeros of their Wronskians. As a result, it remains to bound 
the number of zeros of Wronskians of polynomials of the form IIj=i f'j 1 ^ ■ 
Such a Wronskian has few zeros thanks to a nice factorization property: after 
factoring out some high powers, we are left with a determinant whose entries 
are low-degree polynomials (or sparse polynomials, depending on the model 
under consideration). It is then straightforward to bound the number of real 
roots of this determinant. 



3.1 Derivative of a power 

We use ultimately vanishing sequences of integer numbers, i.e., infinite se- 
quences of integers which have only finitely many nonzero elements. We 
denote this set N^). For any positive integer p, let S? p = {(s\, S2, ■ ■ ■) G 

oo 

(N) | ^i Si = (go f or each p, this set is finite). Then if s is in we 

8=1 

observe that for all i > p + 1, we have Si = 0. Moreover for any p and any 

oo 

s = (si,S2,...) € N^), we will denote \s\ = Yl s i (* ne sum makes sense 

i=l 

because it is finite). 
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Lemma 10. Let p be a positive integer and a > p be a real number. Then 

where (f3 as ) are some constants. 

We define the total order of differentiation of a differential polynomial of a 
function with an example: if / is a function, the total order of differentiation 
of f (f'f (/ (4) ) 3 + 3//' is max(3 *0 + 2*l + 3* 4, 0*1 + 1*1) = 14. 

Lemma [TU] just means that the p-th derivative of a power a of a func- 
tion / is a linear combination of terms such that each term is a product of 
derivatives of / of total degree a and of total order of differentiation p. This 
lemma is proved in Appendix [A] 

3.2 Several models 

In |10j . the authors gave an £°( m2fc ) bound on the number of distinct real 

k m 

roots of polynomials of the form / = X^ilT/i ' ' 3 > w here the fj are poly- 
pi j l 

nomials with at most t monomials. We improve their result in Theorem [12] 
using our results on the Wronskian from Section [2] 

Lemma 11. Let M be a set of T monomials and fi, ■ ■ ■ , f s be polynomi- 
als whose monomials are in M. For every formal monomial P in the s 2 
variables fi, /{, . . . /* ~ , f%, fy, ■ ■ ■ > /J* -1 of degree d and of order of differen- 
tiation e, the number of monomials in x of P(fi, /{,..., is bounded 
by (^y^ 1 )- More precisely, the set of these monomials is included in a set 
Ed,e °f s "i ze a t rnost (^y^ 1 ) which does not depend on P. 

Proof. Let M d be the set of monomials which are the product of d not 
necessarily distinct monomials of M. The cardinal of this set is bounded 
by the cardinal of the set of multisets of size d of elements in M , that is 
(^^T^ 1 )- ^ is easy to see that we can take the set Ed >e defined as the set of 
monomials of x~ e M d . Its cardinal is bounded by the cardinality of M d . □ 

k m 

Theorem 12. Let f = U f ' J be a non identically zero function such 
that fj is a polynomial with at most t monomials and such that m € R and 

ay G N. Then, Z R (f) < Aktm + 4 (e(l + t))^ =0(rf). 

Moreover, if I is a real interval such that for all j, fj(L) C ]0,+oo[ 
(which ensures f is defined on I), then the result is still true for real (possibly 

negative) powers ay, i.e., Zj(f) < Aktm + 4 (e(l + t)) 2 . 



/w*- |si n (/ (fc) ) 
fc=i 
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Proof. Let N an integer such that for all i and j, we have a%j + N > 0. 
Let us consider / = Yli=i a i9i where gi = YYJLi f^ l ' 3+N+k . Note that / = 
/ " Ilj=i fj' We are going to bound the number of zeros of f. In both 
cases (whether otij are integer or real numbers), the functions gi are analytic 
in /. Furthermore, we can assume without loss of generality that the family 
(gi) is linearly independent. Indeed, if it is not the case, we can consider 
a basis of the family (gi) and write / in this basis. Then we can suppose 
that all a, are non-zero, otherwise, we remove these terms from the sum. We 
want to bound the number of zeros of W(gi, . . . , g s ) for all s < k to conclude 
with Theorem [9] We know that for 1 < u, v < s 



g(v-l) 



E n(/; 



ri,r2,...,r m j=l 
ri + ...+r m =v— 1 



'j 



to) 



(4) 



We use now Lemma 10 and we simplify the notation by writing /3 u j s instead 

°f Pa u i+N+k,s- 



si- 1 ' 



e n 

r 1: r 2 ,...,r m j=l 

ri + ... + rm=v — l 



a u>j +N+k-\s\ 



nil 



k=i 



c (fc) 



m 



a uJ +N 




l<p,g<s 



with : 



l<p,q<s 



e n 

ri,ra,...,r m j=l 

ri+...+r m =v— 1 



(5) 



E a,-,./;- 1 -'* 1 n (/ 



fe=l 



The polynomial T n) „ is homogeneous of total degree (t> — l)m with respect 
to the s 2 variables ( /i 9 ^ ) and each of its terms is of differentiation 

order v — 1. 

Then, we notice that, in the first parenthesis does not depend on v 
and the second one on u. We get 



W(g 1: ...,g s 



nn/ 

Vi=ii=i 



Hence, 



dij+N+k-i+l 
j 



det ( (T u , v l(f^-%< p , q < s 



U,V<1 



Z(W( 9l ,...,g s )) < [ ] +Z(det(r u ,,, 



^)l<P.9<s 



• (6) 
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We are now going to bound the number of monomials in x of det(T U]V ). 
We saw that T UjV is a homogeneous polynomial of degree (v — l)m with re- 
spect to the s 2 variables ( fi 9 ^ ) and of order of differentiation v — 1. 

^ > l<P,9<s 

Moreover, as the family (<7j) is linearly independent and as these functions 
are analytic, the Wronskian is not identically zero (Lemma|5|. So det(T u v ) 

is a linear combination, with respect to the variables (fp q ^) , of 



l<p,q<s 
s 

monomials of degree exactly ^ (v — l)m = m,(^) and of order of differentia- 



v=l 



tion Q). By Lemma 11 the monomials in x of each term of det(T UjV ) are in 
the set / s \. Consequently, the number of monomials in x of det(T u v ) 

is bounded by the cardinal of -^( s ) m ( a y i- e - i™ 1 ^ 2 }^-! *)' Descartes' rule 
of signs (Lemma [2]) gives 

We have now all the tools to prove the theorem. We have: 

k 

Z(f) < Z^am). 

i=l 

By Theorem |9l 



Z(f)<k-l + 2j2z(W( gi ,...,g s )). 



s=l 



Using formula ([6]): 



z(f) < k - 1 + 2k [ z (fj) + 2 E z ( W 9 ~ 1] )p.9<' 



By Descartes' rule, YlJLi^ifj) < (2£ — l)m. We can then apply (|7| to 
obtain the inequality 

Z(/) < fc - 1 + 2fc(2t - l)m + 2 E (2 ("^^ ^ " " X ) • 
Finally, we use the well known bound: (^) < (en/k) k 

Z(f) < k - 1 + - 2/cm - 2/c + 4 + 4 V f e [ 1 + mt ~ 



< 4fetm + 4(e(l + *)) 2 . 



□ 
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Using polynomials of small degrees instead of sparse polynomials, the 
same argument gives a polynomial bound. 

Theorem 13. Let f = Yli=i a, i XYj=i fj^'* w h ere / ^ s n °t nu ll, the fj are of 
degrees bounded by d and such that a% are reals and a% are integers. Then, 
Z R (f) < \k 3 md + 2kmd + k ~ ^ 

Moreover, if I is a real interval such that for all j, fj(I) C R + * (which 
ensures f is defined on I), then the result is always true for real powers on, 
i.e. Zi(f) < \k*md + 2kmd + k ~ 



Proof. In the proof of Theorem 12 we saw that det (T uv ) is a homogeneous 

polynomial of degree m{^j in the s 2 variables fx, /(,••• , fs ■ So, it is of 
degree mdfy in the variable x. Moreover, as the family (gi) is linearly inde- 
pendent and as these functions are analytic, the Wronskian is not identically 
zero (Lemma [5]). In Equation Q, the first term is bounded by md and the 
second one by md(^j. By Theorem [sfl^j the number of zeros of 5^f=i a i9i ^ s 
bounded by (^k^md + 2kmd + k) . □ 

Avendaho studied the case / = J^Li x ai (ax + b)P* where cti and are 
integers [1]. He found an upper bound linear in k for the number of roots. 
But he showed also that his bound is false in the case of real powers. We 
find here a polynomial bound but which works also for real powers. 

Corollary 14. Let f = Yli=x CiX ai (ax + b)^ . Let I be the interval {x 6 
R\x > A ax + b > 0}. Then Z^f) = 0(k 3 ). 

Li, Rojas and Wang [16J (Lemma 2) showed that polynomials / = 
Yli=i a i Iljli( c j^ + dj) ai ' J where coefficients cij, bj, Cj and exponents atij 
are real, have at more m + m 2 + . . . + m k ~ l zeros. Our result improves this 
bound: 

Corollary 15. Let f = Yli=i a i YYj=i( c j x + dj) ai - j where the coefficients 
a%,Cj, dj and the exponents ctij are real numbers. On the interval L = {x G 
R|Vj, cjx + dj > 0} we have z\{f) = 0(mk 3 ). 



4 Some Algorithms for Polynomial Identity Testing 

A PIT algorithm takes a polynomial as input, and decides whether the poly- 
nomial is identically equal to zero. There are two classical forms for these 
algorithms: blackbox and whitebox. For the first one, the input is given by a 
blackbox. And in the second case, the input is given by a circuit. These two 
types of algorithms are not comparable in our case since, if we have a circuit, 
we cannot always evaluate it efficiently on an input because the circuit may 
be of high degree. 

1 Using Theorem [8] instead of Theorem [9] would multiply our upper bound by a O(k) 
factor. 
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4.1 Blackbox PIT algorithms 

The bounds on real roots of Theorem [12] immediately give a blackbox PIT 
algorithm for some families of polynomials. 

k m 

Corollary 16. Let f = J^a^ Y\ fj *' J be a function such that fj is a polyno- 

<=l j=i 

mial tmi/i at most i monomials and such that a, € IR and a$j G N. Then, 

there is a blackbox PIT algorithm which makes only 1+Aktm+A (e(l + £)) 2 
queries. 

Proof. We consider the algorithm which tests if the polynomial outputs zero 



12 



this set is 



on the 1 + 4ktm + 4 (e(l + i)) 2 first integers. By Theorem 
a hitting set, that is, if the polynomial is not zero, then at least one of these 
integers will not be a root of the polynomial. □ 

Corollary 17. Let f = Y2i=i a i TljLi /T*'* where the fj are of degree 
bounded by d and such that aj are reals and oti are integers. Then, there is 
a blackbox PIT algorithm which makes only 1 + ^k 3 md + 2kmd + k queries. 



Proof. We apply Theorem 13 □ 



4.2 A Whitebox PIT algorithm 

Results for whiteboxes are more complicated. They rely on the link between 
Wronskians and linear independence. In this section, we prove the following 
proposition: 

k m 

Proposition 18. Let f = Y2 a iTlfj w here fj is a polynomial with at 

i=l j=l 

most t monomials, the are integers and the aij are non-negative integers. 
Let C be an upper bound on the degrees of the fj, on the bit size of their 
coefficients as well as on the bit size of the coefficients a,{ and of the exponents 
cti,j . Then, there is a whitebox PIT algorithm which decides if f is zero in 

time 6 ^C , 2 4mfc2lo § t ). 

First, we will need an algorithm for testing if some Wronskians are 
identically zero or not. We now describe an algorithm which takes as in- 
puts functions h\ = Y\J=i(fj) ai ' j ^ ■■■ ,h>i = Y\jLi(fj) ai,:i (given by sequences 
(fj)i<j<m an d (aij)i<i<! i<j<m) an d which outputs the leading coefficient 
of the Wronskian W(hi, . . . ,hi) if this determinant is not identically zero, 
and outputs zero otherwise. We will use the notation fin) = 0(g(n)). It is 
a shorthand for f(n) = O (g (n) log k g(n)) for some constant k. 

Proposit ion 19. There is an algorithm which on the input (fj)j<m> 
(cti,j)i<l,j<m outputs the leading coefficient of the Wronskian of 



10 



Ibf/r- ii^r- 

v i=i 3=1 
i/ i/ie Wronskian is not identically zero and which outputs zero otherwise. 
This algorithm runs in time O ( C2 Aml2 log * 



Proof. As in the proof of Theorem [12j we in fact compute the Wronskian of 

I™ 



(<7i, . . . , <#) where gi = YYjLi (fj) ■ To get the correct leading coefficient, 
we can just notice that: 



W(g l ,...,g l ) = W(h 1 ,...,h l )ll(f j f . 

3=1 

Hence, we want to compute the Wronskian of (g%, . . . , g{). Again as in 
the proof of Theorem jl2j we factorize each column u by Iljli (fj) au ' j an d 
each row v by Ilj=i {fj) We will denote the resulting matrix by M. 

The entries of this matrix are polynomials. 

According to Lemma [28] in Appendix [B] we can compute the expanded 
polynomial in one cell (v,u) of M in time O (2 vm t mv v m Clogl). 

Then, computing all entries of the matrix which is of size (I x I) needs 
O {2 lm t ml l m C^ operations. Next, we have to compute the determinant of this 
matrix. We are going to compute this determinant directly by expanding 
it as a sum of /! products. This computation takes time O ( C2 4m ' 2 log ' ) by 



Lemma 29 in Appendix [B] 

If the determinant is zero, it means that the Wronskian is zero, then 
the algorithm outputs zero. Otherwise, for computing the leading coeffi- 
cient, we have to multiply the coefficient we got by the leading coefficient of 

(nU nr=i(/j) a ^)(nU wlmY^) This operation takes q ( Cm i(c + 1)) 

operations since we can compute the product of n integers of size s in time 
0(ns). □ 

Second, we will also need the following algorithm: if W(h\, . . . , hi) ^ 
and W(hi, . . . , hi+i) = then find a\, . . . , a;+i such that a\h\ + . . . + aihi = 
hi + \ (these constants exist according to Lemma [5]). So for each i G IN, 
a±h^ + . . . + aihi = hfK. We can compute the a,j using Cramer's formula. 
As a result, for each 1 < j < I we have: 
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hi 
K 



hj-i 



hi- 



l+l 



(i-i) 



h 



(i-i) i,0-i) 



1+1 



n j+ i 



n j+i 



hi 
h' 



lc(W(hi,... 



W(hi,...,hi) 
• hj-i, hj+i, hj+i, . . . , h;)) 



ic(W(hi,...,h,)) 

where lc(W(h\, . . . , hi)) is the leading coefficient of the Wronskian for the 
family (hi, . . . , hi). The previous algorithm (Proposition 19) computes these 
coefficients, so we can compute the (ay) in time O ^C2 4m ' 2 lo s*^ . 

Finally, with such algorithms, we just need to go from a\h\ to a\h\ + . . .+ 
Ofchfc. Each time we add a hj, either it is linearly independent and we add it 
to the current basis or it is dependent, and we write it in the current basis. 
At the end, a\h\ + . . . + a^hk is expressed as a linear combination of basis 
functions. We just have to check if all coefficients are zero to conclude if this 



function is identically zero or not. This completes the proof of Proposition 18 



5 An improved upper bound 

In this section, we give a proof of Theorem [9] First, we point out that 
Voorhoeve and van der Poorten's paper [25j contains a result similar to 
Theorem [8j except that all zeros are counted with multiplicity. 

Theorem 20. [Voorhoeve and van der Poorten, 1975] Let fi, ■ ■ ■ ,fk be real 
analytic functions over an interval I. Then, 

fc-2 fc 
N(h + ... + f k )<k-l + ^N(W j ) + Y J N{W j ) 

3=1 3=1 

where N(f) is the number of roots of f on I counted with multiplicities. 

This result immediately implies Theorem [7] in the analytic case. However 
in our applications we will have to not consider the multiplicity of roots. 
Using some ideas of the proof of this theorem, we can nevertheless improve 
equation 

We will denote W{ = W(fi, /bj • • ■ ; fi) for i > 1 when the family 
(fi, . . . , fi) is clear from the context. Finally, we define Wq = 1. 

In addition to Lemma |5j several connections between the Wronskian and 
the linear combination of the functions are known. We will use a result of 
Frobenius ETJ [TT] : 
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Lemma 21. Let /, a family of analytic functions. Let Ri be the family of 
functions defined by: 

Ro = h + • • • + fk 
W? +1 ( Ri 



Then the functions Ri are analytic and Rk—i = Wk- 

We can now give a proof of Theorem [9] Let Ri be the family of rational 
functions defined by: 

Ro = fx + .-- + f k 



R 



i+1 



Wf +1 ( Ri 



Wi \w i+1 . 

We will prove by induction that for all 1 < i < k — 1, the analytic function 
Ri has at least Z(fi + ... + f k )-i- Z{Wi) - 2 Z(Wj) roots on I. That 

yields the theorem with i = k — 1 and Lemma [21] 

If i = 0, then Rq = fi + ■ ■ ■ + fk and i?o has exactly (and so at least) 
Z(fi + ■■■ + fk) zeros. 

Suppose that the property is verified for a particular i such that i < k — 2. 
We will denote m x {F) the multiplicity of the root x in F for all x G R and 
analytic function i* 7 . (If x is a zero of order fc, we will denote m x {F) = — k). 
We define four values (roots are counted without multiplicity): 

• Z^~ is the number of x G R such that m x ( i ) > ^^(Wi+i) > 0. 

• is the number of x E R such that m i; ( ^ i ) = m x (Wi+i) > 0. 

• ^ r> is the number of x G R such that m x (Wj+i) > m^ y^ ) > 0. 

• is the number of x G R such that m^IVj+i) > = m x( 

We have: = Z r + Z = + Z+° + Zf>. 

We know by Lemma 21 that is indeed analytic. Then by induction 
hypothesis, the rational function ^ has at least 

Z(fi + . . . + f k ) - i - Z(Wi) - 2 HWj^J - Z= - Z7> 

roots and at most Z^ y + Z~° poles. By Rolle's Theorem, the number of 
zeros of is at least 



■1-1 



Z(fi + ... + f k )-i- Z(Wi) - 2 [ ^2z(Wj) )-Z=- Z s 

3=1 



-> 



[Zr> + Zf + 1] 
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So, the number of zeros of R-i+i = ^ ( ™^ 1 at least 



-> 



Ztfi + ..- + fk)-i- Z(Wi) - 2 \J2 Z(Wj)j - Z= - Z t 

-[Zr> + Z r° + i} + zr>-Z(W t ). 
We used here that if x is such that < m x {Ri) < m x (Wi+i) then 



-m x (W i+l ) < m x < 



and so m T 



w * + i f Ri V 



> m x (W i+ i) 



1 > 0. Hence 



Z(R i+l ) > Z(fi + . . . + f k ) - (i + 1) - Z(W i+1 ) - 2 i J2 Z(Wj 
6 Optimality of Theorem [8] 

Recall that in Theorem |8j it was proved that 

Z(ai/i + ... + o fc / fc ) < (l + |T|)fc 
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that this 



with T = |J Z (W(/i, . . . , fi)). It will be shown in Theorem 

l<i<fc 

theorem is quite optimal in the sense that for arbitrarily large values of |T| 
and k, we can find functions fi, ■ ■ ■ ,fk an d coefficients a±, . . . , such that 

Z(oi/ a + . . . + a fc / fc ) > (l + |T|)(fc-l). 

We begin by proving a technical lemma. 

Lemma 22. Let f be a non-constant polynomial. There exists, for < 
i < q, rational functions F^ q such that if we define the function h Pt i = 

(^JT (f'Y P~\ then we have: 



1. for all q > 0, we have F t 



1 



2. for all < % < q, the rational function [(f) Fi,q\ is a polynomial 



3. For all q > 0, for all p > 1 we have: (f p ) 



(<?) 



i=0 



The main point is that F^ q does not depend on p. 

Proof. We define Fi A by induction on q. If q = 0, let us define i^o = 1- 
Then, we have f p = h Pj o and (/')° Fq 
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We suppose now that F itQ i are denned for all i, q' such that < i < q' < q. 
Let us define We have: 



(h p ,i)' 



(p-i)\ 



^ [(? - iKfT 1 ! 13 - 1 - 1 + ^f"(fy- 1 f p - ^ ] 



f.f"\ 

- h p , i+ i + I i— \ h Pji . 



So, 



(n (q+1) = [E h p^ 



^^=0 



-E 



hp,i(F ijq )' + ( h Pi i+i + ( ) /i P ,i ) 



h p ,o{Fo,q)' + 



i-l,q 
+ hp tq+ lF qtq . 

We can then define 

for 1 < i < q, F iiq+ i = (Fi, q )' + F i>q + Fi_i, q 

and Fg+i^+i = = 1. 

Then, we have (1) and (3) by construction. Finally, (2) is verified, since 
by induction hypothesis: 

U') q+1 K q = [/' {(f') q Fi, q )]' - (? + 1)/' ((/T^s) ^ a polynomial. 

□ 

We are going to show that the zeros of W (f ai+k , . ■ ■ , f as+k {x)) are either 
zeros of / or zeros of /'. 

Lemma 23. Let f be an analytic non- constant function in an interval I and 
ao, . . . ,«fc be k pairwise distinct integers (with k>\). Then 

[x eI\3s<k,W (r i+fc , . . . , f a ° +k ) (x) = 0} C {x G I\(ff')(x) = 0} . 
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Proof. Let us consider f ai+k , . . . , f ak+k . First suppose that this fam- 
ily is linearly dependent. This means that there exist some constants 
(oi,...,a fc ) G R fc \{(0,0, ...,0)} such that 



^o 4 f' +fc = 0on I. 



(8) 



i=l 



But the integers {oci + k) are all distinct so the polynomial P(Y) = 
Ei=l is not zero. Hence P(Y) has a finite number of roots. By 

Q, Im(f) is included in the (finite) set of roots of P. Nevertheless, as f is 
continuous, by the intermediate value theorem Im(f) is a real interval. So 
Im(f) is a singleton. This contradicts the hypothesis that / is not constant. 
Therefore, this family is linearly independent. 

Let A be the matrix defined by Ay = t^" 1 ) 
J-1 

get Ajj = ^2h a . + ki F[j_i, i.e. in terms of matrix product: 



. By Lemma 
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we 



A - [^oi+fc,I-l]i<i,i< a [ F «-iJ-i%i]i</,j<s • 

The second matrix of the product is an upper triangular matrix whose entries 
on the main diagonal are 1 and so its determinant is 1. Then, 



det l(Ay) 1 



<i,j<s 



det ((^o i +fcj-i) 1 < iii <,) 



Finally, h ai+k>j ^ = [/*] [(/' 1 

does not depend on j and the second one on i. Consequently, 



The first bracket 



J2 a l 



(/')(a) /( fc - a + 1 >+(£) 



det (h ai+kt j-i) 
Then, for all x in I : 

W ( f ai+k f&s + l 

f(x) = or /'(x) = or det 



det 



(ai + fc)! 



(ai + A;- j + l)7 ' 



^» det (/ia i+fe j-i) (a) = (9) 

(q, + fc)! 
(oi + fc- j + 1)! 







If det 



(ai+k)\ 



0, as it does not depend on x, the function 



(ai+k-j+iy., 

det (h ai+ kj-i) vanishes for all x and so the Wronskian is zero over /. But 
as the functions f ai+k are analytic, they would be linearly dependent by 
Lemma [5] That contradicts the hypothesis. Consequently, 

[x G I\3s <k,W . . . , r +k ) = o] C {x G = 0} . 



□ 
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As a byproduct, wegive another proof of the weak version of Descartes' 



rule of signs (Lemma 
number of distinct rea 
with f(x) - 



Lemma 
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21) . Let g = X/i=i a i x T ■ We nee d to show that the 
roots is bounded by 2k — 1. We can use the result of 



x. In this case g = Y^,i=i a iX a ' '■ We get 
T = {x G I\3s <k,W . . . , n +k ) = {0}} C {x G I\(ff')(x) = 0} 

So, Theorem [8] gives: Z(f) < 2k — 1. A similar proof of Lemma [T] appears 
in [20] (Part V, exercise 90). 

In Lemma [23] it can be seen that the converse of the implication (J9| is 

true as soon as /' really appears as a factor of det ^(/i ai+ ^j_ i) 1<i j < A ■ It 

is the case when (X) is different from zero, that is to say when s > 2. This 
implies the following result. 

Lemma 24. Let f be an analytic non-constant function in an interval I and 
ato, . . . ,ctk be k pairwise distinct integers with the condition k > 2. Then 

[x G I\3s <k,W (r i+fc , . . . , f a ° +k ) (x) = O} = {x G I\(ff)(x) = 0} . 

We have now all the tools to prove the main result of the section: the 
optimality of Theorem [8] 

Theorem 25. Let T = {x G L\3i < k, W(/i, . . . , fi){x) = 0} as in Theorem 
[s| For every k and p, there exists a function g = X^=i a if ai such that a, 
are positive integers, f is a polynomial such that |T| > p and such that g has 
at least (1 + |T|) (& — !) + Z(f) zeros. 



) . This polynomial is /c-sparse and has 2k — 1 
. < -1< < 1< k- 1. 



Proof. Let h = x Y\ k = i (x 2 — 
distinct real roots: — k + 1 < 

Let / = fc ]1 (x ~ 2t). 

i=l 

Then, we just have to verify that g = ho f has the required properties. 
We have g(x) = if and only if f{x) G [—k + 1, /c — 1] n Z. But for y an 
odd integer, we have \ f(y)\ > k — 1 and for y an even integer between 2 and 



2 + 2 



2+1 

2 



we have f(y) = 0. By the Intermediate Value Theorem, g has 



at least k — 1 zeros over each interval (n, n + 1) with 1 < n < 2 + 2 
So 



p+1 

2 



2 11 + 

(2fc 



1 



"p+1 
2 

1 + 



(A; — 1) + ( 1 + 



(10) 
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Rolle's Theorem ensures that for two roots of /, there exists a root of 
/' which is strictly between both roots of /. Hence, Z(ff') > 2Z(f) — 1 = 



1 + 2 



p+1 

2 



. Considering the degree of //', we find Z(ff') = 1 + 2 



p+i 



Besides, / is not constant so by Lemma 24 |T| = Z(ff'). Hence, 
|T| = Z(ff') = 1 + 2 P+ 1 
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We can verify that the hypothesis |T| > p is true. Finally, equations (10) 



and (11) show that Z(g) > (|T| + l)(Jfc - 1) + Z(f). 



□ 



In the proof of Theorem 25 the roots of all W(f\, . . . , fi) are included in 
the zeros of W(/i, . . . , /&). So, it could be possible to improve both Theorem 
[8] and Theorem [9] by proving the following proposition. 

Open question. Let fx, ■ ■ ■ , fk be analytic functions on an infinite interval 
I and ai, . . . ,Ofe be non-zero real constants. Is the inequality 

k 

Z(oi/i + • • • + a k f k ) < k - 1 + £ Z(W(/ la . . . , /i)) 

i=l 



always true? 
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A Proof of Lemma 10 



In the following, is the sequence (0,0,..., 0, 1, 0, 0, . . .) where the 1 appears at 
the i th coordinate. We show this lemma by induction over p. If p = 1, then 
(/<*)' = af f a - x . That is the basis case since S*i = {(1,0,0, . . .)}. We notice that 

Pa, (1,0,...) = Ot. 

Let us suppose that the lemma is true for a fixed p. By induction hypothesis, 
we have 




£ ^-r ' IK/ (fc) )" 



k=l 
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where 



. 91 =EM^ N )'(n (/ ( T) 
32= e ^f a ~ is{ { n (/'• 

By rewriting each term, we get 

.91= E ^(a-ki)/'/^ 151 - 1 ]!^) 8 ' 
= e /3^(«-k'i + i)r- |s,| n(/ (fe) 



fc=l 



.92 



E/3a, s r^E^ +1) ' 1 1 (-/ 



(AO 



E 



p+1 



E /W ej -e J + J aHS,| (4 + l)Il(/ 



s-e,.+e 



p(fe) 



k=l 



If s is in S^ p , then s+ei € ^p+i and if moreover Sj ^ then s— ej+ej + i € =5^ p +i. 
So the result is proved and the constants f3 are defined by: /3 a ,(1,0,0,...) = a ; an d if 
s € S p with p > 1, then 

/3a, s = l Sl ^o(a - |«| + l)/?a,(si-l,S2,S3,...) 

+ E (^i- 1 + 1 )/3c t ,(si,...,s 3 _ 1 ,s J _i + l,s J -l,s j + i,...)- 
j ■ 2<3<P 

B Proof of lemmas for Section 0] 



In this section we prove two lemmas ( |28| and [29| needed for the proof of Proposi- 
tion [19] These proofs are elementary but somewhat technical. 

In the following, we will compute additions of n integers of size s in time 
0(n(s + logn)) and products of n integers of size s in time <D(ns). 

We begin by bounding the complexity of expanding a product of sparse poly- 
nomials. 

Lemma 26. We consider a product of /i r '-sparse polynomials Pi, . . . , P^ of degrees 
bounded by 7 with integer coefficients of size bounded by 7. This product can be 
expanded in time 0(r At 7). Moreover, the size of new coefficients is bounded by 
^7 + ^logr. 



21 



Proof. For expanding such a product, we compute one by one each monomial of 
the sum and we store the coefficients of these r M new monomials. For computing 
one coefficient, we have three things to do. We have to compute its degree (sum of 
/i integers of size 7) in time O (fj,( / y + log /i)), its coefficient (product of fj, integers of 
size 7) in time O (fJ-'y) and we add together the monomials with the same exponent. 
At the end, at most r M coefficients will be added together to form a given monomial, 
so the size of the coefficient is bounded by nj+fi log r. Hence, as we add coefficients 
one by one, at each step we have to add an integer of size prj by one of size at most 
/17 + filogr. Each term of the sum takes time 

6 (n(-y + log fi) + fij + (M7 + ju log r)) 

= 6 + ^log(flT)) . 

Therefore, computing all coefficients takes time 

6(T"(A*7 + Mlog(MT))) 

= 6 (7^7) . 

□ 



Theorem 12 uses some constants (3 a s which have been defined in Lemma 10 



We will need to compute them. 

Lemma 27. For every p in N, we have \S p \ < 2 P ~~ 1 . For every a, p in N and for 
every I in S p , < f3 a s < (p 2 + a) p . 

Furthermore, for every a, p in N we can compute all /3 ajS with s € S q and q < p 
in time 0(2 p loga). 



Proof. We showed in the proof of Lemma 10 that (3 a ,(i,o,o,...) = a an d if s G S p 
with p ^ 1, then 

Pa,s = lsi^o(a - |«| + 1)^0,(81-1,82,83,...) 

+ ^2 ( S J-1 + 1 )/ 3 a,(si,...,8 i -i,s J _i+l,s 3 — l.sj+i,...)- (12) 

3 ■ -'• r ;■ 

Sj7 t0 

However, in the formula above, the sequences (si — 1,S2,S3,...) and 
(si, . . . , Sj_i, Sj_! + 1, Sj — 1, Sj +1 , . . .) fall in S p ^i. Let us denote M a p = maxl/Ja sl. 

Hence, if p ^ 1, (12 1 implies, 

M ayP < (p 2 + a)M QiP _i. 

Since, M a .i — a, we get by induction /3 ajS < (p 2 + a) p ~ 1 a. 
For computing these constants, we notice that: 

- \{s g 5 p+1 | Sl ^ 0}| + |{a G 5 p+1 | Sl = 0}| 
<2|5 P |. 

The inequality comes from the two surjective functions: 

S p —> {s G S p +i\si ^ 0} and S p —> {s G 5 p +i|si = 0} 
s H' (si + l,s 2 , . . .) S4 (0, s 2 , ■ • ■ ,s Sl ,si +Sl + l,s 2+Sl , . . .). 
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Hence by induction, |«S P | < 2^ and \{} v S p \ = 0(2 P ). Then, if s € S p is fixed, 



for computing f3 a s with (12 1, we need to compute p products of an integer of size 
plog(p 2 +a) by an integer of size logpor log a (in time O (plog(p 2 + a))) and asum 
of all these products in time (p 2 \og(p 2 +«))■ Finally computing all constants 
Pa,s with s 6 S q and q < p needs time 0(2 p p 2 log(p 2 + a)). That proves the 
lemma. □ 

We can now prove the two intermediate lemmas of Proposition [19] For the 
following, we keep the notations of Proposition |19| 

Lemma 28. Computing the expanded polynomial of a cell (v,u) in the matrix M 
takes time O (2 vm t mv v m C\ogl). Coefficients are of size bounded by O (mvC log tl) 
and degrees of size bounded by Cmv. 



Proof. We also keep the notations of Theorem 12 Each cell (v,u) corresponds to 
the polynomial 



T, 



in 



l<p,q<l 



E 

r-i,ri,...,r m j = \ 

n+...+r m =»-i 



n 



E 



k=l 



E E 



r\,ra,...,r m s 1 ,s 2 

1 +...+r- m =v-l 



II v. II /; - 



c(fc) 



\3 = 1 



k=l 



The first sum is a size at most v rn and the second one of size bounded by 2 vm 
(Lemma 27 1. For computing one term, first, we need to compute Yij=i0u,j,si — 
Ilj=i @a u j+i,si wich is a product of m integers, each one of size v\og(v 2 + 2 C + I) 
(since a < 2 C ). It is done in time 0(mvC\ogl). 

Now, we want to develop the formula with respect to x. We saw in the proof of 

that each monomial with respect to the I 2 variables (jp q 1 ^ 



Theorem 



12 



l<p,q<l 



of total degree m(v—l). We consider one monomial with respect to (fp Q 1 M 

V / l<p,q<l 

It is a product of m(v — 1) t-sparse polynomials with respect to the variable x. By 
Lemma 26 this product (the second parenthesis in the formula) can be expanded 
in time O (f m (" -1 ) C*) and coefficients are of size Cm(v — 1) + m(v — 1) logi. Then 
each coefficient is first, multiplied by the corresponding coefficient 0.7=1 &u,j,si m 
time 

6 (max{Cm(u - 1) + m(v - 1) log t, mv \og(v 2 + 2° + I)}) 
= 6(mvC\ogtl) , 

that gives an integer of size at most O (mvClogtl). Second, it is added to the 
stored coefficient corresponding to the same monomial in time 

6 (mvClogtl + log(v m 2" m )) 
= O (mvClogtl) . 
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To conclude, computing the cell takes time 

6 fv m 2 vm ^ mvC \ogl + (t m(v -^C + t m(v - 1 \mvC\ogtl + mvC\ogtl) 
= O (2 vm v m t m(v -VC log ifj . 

That completes the proof of the lemma. □ 

Lemma 29. Assume that the entries of M are given in expanded form (i.e., as 
sums of monomials) . Computing the determinant of M takes time O (C2 4ml 
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Proof. For each one of the (7!) permutations, we expand the corresponding poly- 
nomial. Each cell has at most 2 ml l m t ml monomials. Hence, each permutation 
corresponds to a product of size I of (2 m '? ?n <" l ')-sparse polynomials. Powers are 
bounded by Cml and coefficient sizes are bounded by 0(mlC log t). By Lemma 
each permutation can be computed in time 

6(2 ml2 r l t ml2 mlClogt) 
= d( 2 3ml2 



and size of coefficients is bounded by 

6 (mPClog t + llog (2 ml l m t ml )) 
= 6 (mfClogt) . 

For computing the whole determinant, we compute permutations one-by-one, 
adding each time new coefficients to the one computed before. This is done in time 



O (l l [2 3ml i °e t C + 2 ml r l t ml ml 2 C\ogt^ 

□ 
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